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We propose a multiterminal nanostructure for electrical probing of the quantum spin Hall effect 
(QSHE) in two-dimensional (2D) topological insulators. The device consists of a ferromagnetic 
(FM) island with precessing magnetization that pumps (in the absence of any bias voltage) pure 
spin current symmetrically into the left and right adjacent 2D TIs modeled as graphene nanoribbons 
with the intrinsic spin-orbit (SO) coupling. In the reference frame rotating with magnetization, the 
device is mapped onto a DC circuit with twice as many terminals whose effectively half-metallic 
ferromagnetic electrodes are biased by the frequency huj/e of the microwave radiation driving the 
magnetization precession at the ferromagnetic resonance conditions. The QSH regime of the six- 
terminal TI|FM|TI nanodevice, attached to two longitudinal and four transverse normal metal 
electrodes, is characterized by the SO-coupling-induced energy gap, chiral spin-filtered edge states 
within finite length TI regions, and quantized spin Hall conductance when longitudinal bias voltage 
is applied, despite the presence of the FM island. The same unbiased device, but with precessing 
magnetization of the central FM island, blocks completely pumping of total spin and charge currents 
into the longitudinal electrodes while generating DC transverse charge Hall currents. Although these 
transverse charge currents are not quantized, their induction together with zero longitudinal charge 
current is a unique electrical response of TIs to pumped pure spin current that cannot be mimicked 
by SO-coupled but topologically trivial systems. In the corresponding two-terminal inhomogeneous 
TIjFM|TI nanostructures, we image spatial profiles of local spin and charge currents within TIs 
which illustrate transport confined to chiral spin-filtered edges states while revealing concomitantly 
the existence of interfacial spin and charge currents flowing around TI|FM interfaces and penetrating 
into the bulk of TIs over some short distance. 



PACS numbers: 73.43.-f, 72.25.Pn, 76.50.-l-g, 72.80.Vp 

I. INTRODUCTION 

The recent theoretical predictionsii^i^ for the quantum 
spin Hall effect (QSHE) have attracted considerable at- 
tention by both basic and applied research communities. 
In the conventional SHE,^ which manifests in multiter- 
minal devices^ as pure spin current I^' in the transverse 
electrodes driven by longitudinal unpolarized charge cur- 
rent in the presence of intrinsic (due to band structure) 
or extrinsic (due to impurities) spin-orbit (SO) couplings, 
the spin Hall conductance Gsh = 1^' can acquire any 
value {y is the bias voltage applied between the longitu- 
dinal electrodes) 1^1^ Conversely, Gsh = 2 x e/47r becomes 
quantized in four-terminal devices that exhibit QSHE.^ 

The QSHE introduces an example of a new quan- 
tum state of matter — the so-called topological insulator 
^'jij^j6j7j8 jj-^ ^.^Q dimensions — which is a band insulator 
with a usual energy gap in the bulk, but which also ac- 
commodates gapless spin-polarized quantum states con- 
fined around the sample edges. Unlike closely related 
quantum Hall insulators;^ where the bulk energy gap and 
edge states appear due to an external magnetic field, TIs 
are time-reversal invariant systems whose intrinsic SO 
coupling opens a bulk gap while generating the Kramers 
doublet of edge states. These edge states force electrons 
of opposite spin to flow in opposite directions along the 
edges of the sample. Since time-reversal invariance en- 
sures the crossing of the energy levels of such peculiar 



chiral and spin-filtered (or "helical"—!^) edge states at 
special points in the Brillouin zone, the spectrum of a 
TI cannot be adiabatically deformed into topologically 
trivial insulator without such statesi^!^ 

The recent experimentsi^iii on HgTe quantum wells 
have confirmed some of the anticipated signatures^ of 
QSHE in line with transport taking place through heli- 
cal edge states, such as: (i) reduction of the two-terminal 
charge conductance; (ii) its independence on the sample 
width; and (Hi) sensitivity to an external magnetic field 
that destroys the TI phase. Nevertheless, this is still per- 
ceived as an indirect and incomplete detection because it 
did not confirm that conducting electrons along the edge 
were spin-polarized. The very recent nonlocal transport 
measurements on multiterminal HgTe microstructures in 
the QSH regime suggest that charge transport in these 
devices occurs through extended helical edge channels^ 
since their results can be explained via simple applica- 
tion of the Landauer-Biittiker scattering theory of mul- 
tiprobe conductance to TIs attached to several metallic 
electrodes.-i^ 

Undoubtedly, the most convincing evidence for QSHE 
would be to measure quantized Gsh, as the counterpart 
of the quantized charge Hall conductance^ in the inte- 
ger QHE exhibited by two-dimensional electron gases 
(2DEGs). However, this is quite difficult since pure spin 
currents'''^ typically have to be converted into some elec- 
trical signal to be observedJ^ii^ii^ Thus, the key issue 
for unambiguous QSHE detection, as well as for the very 
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FIG. 1: (Color online) (a) The proposed TI|FM|TI het- 
erostructure consists of a central ferromagnetic island with 
precessing magnetization which pumps pure spin currents to- 
ward the left and the right 2D topological insulators modeled 
as graphene nanoribbons with intrinsic SO coupling. In the 
rotating frame in which magnetization is static, device (a) is 
mapped onto a four-terminal DC circuit in panel (b), whose 
effectively half-metallic ferromagnetic electrodes have electro- 
chemical potential shifted by ±ftaj/2 with respect to the Fermi 
level Ef of unbiased normal metal electrodes in the labora- 
tory frame. 



definition of what constitutes direct experimental mani- 
festation of 2D TIs, is to design^, devices where electrical 
quantities can be measured that are directly related to 
helical edge state transport and the corresponding quan- 
tization of GsH- 

The experimentsi^ii^ii^ on the inverse SHE — where 
injection of pure spin current into a device with SO 
couplings leads to deflection of both spins in the same 
direction and corresponding Hall voltage between lat- 
eral sample boundaries or charge current in the trans- 
verse electrodes — provide guidance for QSHE probing 
via conventional electrical measurements. The devices 
constructed for these experiments are very flexible and 
allow for multifarious convincing tests confirming SHE 
physics. In particular, one of the inverse SHE experi- 
ments has injected pure spin current, generated by spin 
pumping^I from a ferromagnetic (FM) layer with precess- 
ing magnetization driven by RF radiation at the ferro- 
magnetic resonance (FMR) conditions, into a metal with 
SO couplings to observe the transverse Hall voltageJ^ 

Here we propose a spin pumping-based nanostructure, 
illustrated in Fig. [ija), where a FM island with pre- 
cessing magnetization pumps, in the absence of any ap- 
plied bias voltage, pure spin current into two adjacent 
graphene nanoribbons (GNR) with intrinsic SO coupling 
that act as the simplest model of 2D TlJu^ This setup has 
an advantage^'' over other possible sources of pure spin 
currents because it evades the conductivity mismatch** 
between metallic injector and TI that would play a detri- 
mental role when spin injection is driven by a bias volt- 
ageJ^ii^ The three-layer central sample is attached to two 
metallic electrodes, and we also analyze six-terminal se- 
tups where additional four transverse electrodes, two per 



each GNR region, cover a portion of its top and bottom 
edge. 

The nonequilibrium Green function (NEGF) picture 
of pumpingfi^ii^ which was utilizedi^ to explain very re- 
cent experiments^ on spin pumping across a band in- 
sulator within a magnetic tunnel junction, converts the 
complicated time-dependent problem posed by the device 
in Fig. [IJa) into a four-terminal DC circuit in Fig. [T{b) 
in the rotating reference frame. This picture also moti- 
vates the usage of a symmetric TI|FM|TI nanostructure 
since in asymmetric devices the central FM island would 
pumpi^ concomitantly a small charge current into the 
TI regions (with quadratic frequency dependence as op- 
posed to dominant pumped spin current which is linear 
in frequency in the adiabatic limit). Within this frame- 
work we obtain the following principal results: (i) in the 
two-terminal device in Fig.[TJa), pumping generates both 
spin and charge local currents inside the TI, but only to- 
tal spin current is non-zero (Figs. [3] and 3]); (ii) imaging 
of local quantum transport in this two-terminal device 
also demonstrates the existence of interfacial spin and 
charge currents at the TI|FM boundary, which are able 
to penetrate into the bulk of TIs over some short distance 
(Fig. [3]); and (m) in the corresponding six-terminal de- 
vice whose TI regions are brought into the QSH regime, 
pumped total spin and charge currents vanish in the lon- 
gitudinal leads, while non-zero charge currents Ip emerge 
in the four transverse electrodes (p=3-6) as the manifes- 
tations of the inverse QSHE (Fig. [S]). The charge con- 
ductances Gt = elp/huj, however, are not quantized. 

It is worth recalling that 2D TIs, also denoted as QSH 
insulators, were initially studied as isolated infinite ho- 
mogeneous systems."^ On the other hand, experiments re- 
quire to embed such materials into circuits where they 
will be attached to multiple metallic electrodes serving 
as current or voltage probes. For example, the very re- 
cent analysis has predicted highly unusual spin dynamics 
at the Tllnormal-metal^i and TI|FM2i interfaces. This, 
together with our findings of interfacial currents around 
TI|FM junction highhghts the need to understand op- 
eration of inhomogeneous nanostructures, with multiple 
TI I metal interfaces and helical edge states of finite ex- 
tent,^^ as a prerequisite for the design of anticipated 
spintronic devices^^ exploiting TIs. Similar issues, al- 
beit without spin dynamics, were encountered in the well- 
known experiments^'^ on inhomogeneous (e.g., containing 
a potential barrier) multiterminal quantum Hall bridges 
where a simple picture of edge transport becomes insuf- 
ficient and one has to take into account trajectory net- 
work22i2i connecting edge channels across the bulk of the 
device;^ 

The paper is organized as follows. In Sec. |TT1 we dis- 
cuss the effective Hamiltonian of the device and NEGF 
approach to computation of pumped spin and charge 
currents in the picture of the rotating reference frame. 
The images of local pumped spin and charge currents are 
shown in Sec. Illll for unbiased two-terminal device, while 
total quantized spin Hall currents in biased and pumped 



charge currents in the transverse electrodes of unbiased 
six-terminal devices are studied in Sec. |TV1 Section IVl is 
devoted to explaining the origin of non-quantized trans- 
verse charge currents and zero total spin and charge cur- 
rent in the longitudinal electrodes of our inverse QSHE 
device using the picture of possible Feynman paths^^ of 
spin-polarized injected and collected electrons within the 
12-terniinal DC device in the rotating frame. We also 
support the conjectured network of such paths using spa- 
tial profiles of local currents computed for the 12-terminal 
DC device. We conclude in Sec. |Vll 

II. NEGF APPROACH TO SPIN PUMPING IN 
MULTITERMINAL QSH SYSTEMS 

The central TI|FM|TI region of the nanostructure in 
Fig. [T] can be described by the effective single 7r-orbital 
tight-binding Hamiltonian in the laboratory frame: 

- l^^cj + ^7so ^^'^ ■ ("^feJ ^ d,fc)cj-. (1) 

Here c\ = (c\^ , ) is the vector of spin-dependent elec- 
tron creation operators and a = (dx, cry,c^z) is the vector 
of the Pauli matrices. For simplicity, the usual (spin- 
independent) nearest-neighbor hopping is assumed to be 
the same 7^; = 7c = 7gnr — 7 for the square lattice 7s; 
of FM island and the semi-infinite leads, as well as for the 
hopping 7c that couples different regions of the device. 
On the honeycomb lattice of TI regions, 7gnr — 2.7 eV 
can reproduce^^ ab initio computed band structure very 
close to the charge neutral Dirac point {Ep = 0). 

The coupling of itinerant electrons to collective mag- 
netic dynamics is described through the material- 
dependent exchange potential , which is non-zero only 
in the FM island. The on-site potential can accom- 
modate disorder external electric field, or band bot- 
tom alignment. The Hamiltonian ([T]) is time-dependent 
since the spatially uniform unit vector m(t) along the lo- 
cal magnetization direction is precessing steadily around 
the z-axis with a constant cone angle 9 and frequency 
/ = uj/2i:. The central TI|FM|TI region is attached to 
two or six semi-infinite ideal (spin and charge interaction 
free) electrodes, modeled on the same square lattice as 
the FM island, which terminate in macroscopic reservoirs 
held at the same electrochemical potential = Ep. 

The third sum in Eq. ([1]) is non-zero only in the 
GNR regions where it introduces the intrinsic SO cou- 
pling compatible with the symmetries of the honeycomb 
lattice. The SO coupling, which is responsible for the 
band gapi Aso — 6\/37so, acts as spin-dependent next- 
nearest neighbor hopping where i and j are two next- 
nearest neighbor sites, k is the only common nearest 
neighbor of i and j, and d^j. is a vector pointing from 
k to i. 
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FIG. 2: (Color online) Illustration of the definition of bond 
spin and charge currents between nearest neighbor sites of 
the square lattice of FM island and nearest neighbor or next- 
nearest neighbor sites of the honeycomb lattice on which the 
2D TI Hamiltonian is defined. The vertical dashed lines show 
cross sections within FM or the TI region where the bond 
currents passing through them are summed up to get the total 
currents in Fig. |4]or conductances in Fig. |6l 



Among several proposals for QSHEf^iii graphene with 
7SO 7^ stands out as the simplest model employed to in- 
troduce QSHE phenomenology,— Z2 classification of TIs,— 
and manifestations of QSHE in realistic multiterminal de- 
viceSj^^i^ as well as to analyze disorder^''^^'^^ and interac- 
tion^i^^ effects on TIs. However, recent scrutiny finding 
minuscule Aso — 1 /^V through density functional the- 
ory (DFT) calculations^^ would push the requirement 
/cfiT <C Ago for the observation of QSHE in graphene 
toward unrealistically low temperatures^^ T ^ 10^^ K. 
However, these calculations are not fully first-principles 
and all-electron — a direct all-electron (i.e., no pseudopo- 
tential) DFT approach^^^ has estimated Aso - 0.6 K, 
which makes the device in Fig. [T] experimentally real- 
izable. The physical explanation of an apparent con- 
troversy between extremely small intrinsic SO coupling 
found in Refs. [H and non-negligible values computed in 
Refs.'a^'ai' is in the fact that 96% of the intrinsic SO split- 
ting in graphene originates from the usually neglected d 
and higher orbitalsi^ 

We select the following device parameters in the anal- 
ysis and Figures below: 7S0 = O.O37, A = 0.67, 
Ep = 10-^^7, / = 2 GHz and = 10°. The ge- 
ometry and size of zigzag GNR (ZGNR) and FM re- 
gions of the device shown in Fig. [T] is characterized by: 
iVj-ZGNR is composed of Nz zigzag chains so that its 
average width is W = aV5{Nz - l)/2 (a = 2.46 A is 
the lattice spacing of the honeycomb lattice); dcNR is 
the number of atoms along the zigzag chain defining its 
length L = (cZqnr — l)a/2; and dpM is the thickness of 
the FM island. For illustration of N^, dcNR, and dpM 
parameters see Fig. [2] 

Note that the selected valuei for 7S0 — O.O37 is much 
larger than the one that would be fitted to first-principles 
calculations The reason for this choice (in fact, even 
larger values have been employed in recent studies2^i21i2^) 
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is the usage of the effective tight-binding Hamihonian ([T]) 
which in the case of vastly different energy scales 730 ^ 7 
would make quantum transport calculations insensitive 
to the presence of 730- On the other hand, this does 
not affect any conclusions about quantized spin trans- 
port properties of a graphene-based model of a TI since 
they do not depend on the particular value of 730 and 
simply require to perform measurements at temperatures 
ksT ^ Aso where the band gap of the TI is visible while 
its Fermi energy is within such gap. 

Although widely-used scattering theory^^ of adiabatic 
spin pumping by FM| normal- metal (FM|NM) interfaces, 
typically combined with the spin-diffusion equation and 
magnetoelectric circuit theory, cannot handle nanostruc- 
tures containing insulators of band (due to spin accu- 
mulation not being well-defined in them) or topological 
type (due to necessity to take into account details of 
one-dimensional transport through helical edge states), 
the NEGF approach to spin pumping^^^iS, can describe 
both cases by taking the microscopic Hamiltonian ^ 
as an input. The unitary transformation of Eq. ([T]) 
via C/ = e*'^'^^*/^ [for m(<) precessing counterclockwise] 
leads to a time-independent Hamiltonian in the rotating 
frame : ^^1^^ 

Hrot = UH,,kit)U^ - ^W—U^ = 7?iab(0) - —a,. (2) 

The Zeeman term hioa^j^^ which emerges uniformly in 
the sample and the NM electrodes, will spin-split the 
bands of the NM electrodes, thereby providing a rotating 
frame picture of pumping based on the four-terminal DC 
device in Fig. [ijb). This term breaks time- reversal in- 
variance (while conserving spin 6*2), but for typical FMR 
frequencies^*' of the order of / ~ 1 GHz, it is smaller 
than^ai^i Aso - 50 \jN . 

The basic transport quantity for the DC circuit in 
Fig. [Ijb) is the spin-resolved bond charge current^i^^ car- 
rying spin-cr electrons from site i to site j 



(3) 



This is obtained in terms of the lesser Green function'^^ 
in the rotating kamo^^ G<{E). Unlike G<{t, t') in the 
laboratory frame, depends on only one time variable 
T = t — t' (or energy E after the time difference r is 
Fourier transformed^). This yields spin 



2e V'^ 



J, 



(4) 



hopping 7ij ^ 0. The definition of bond currents is illus- 
trated in Fig. [21 both for the square lattice of FM and 
the honeycomb lattice of the chosen 2D TI model. 

The rotating frame four-terminal device in Fig. [T](b) , 
or its twelve-terminal counterpart originating from a six- 
terminal QSH bridge [illustrated by insets in Fig. [5], 
guides us in constructing the NEGF equations for the 
description of currents flowing between their electrodes. 
The electrodes in the rotating frame are labeled by (p, a) 
[p = L,R and a =1, |] and they are biased by the electro- 
chemical potential difference /Lt^ — fi^, = hw. Thus, these 
electrodes behave as effective half-metallic ferromagnets 
which emit or absorb only one spin species. 

The rotating frame retarded Green functioni^ii^ 



&\E)^[E~li,,t-^''{E)]' 



and the lesser Green function 



G<{E) ^ G''{E)t.<{E)GHE), 



(6) 



(7) 



describe the density of available quantum states and how 
electrons occupy those states, respectively. Here Hrot is 
the matrix representation of Hrot © in the basis of lo- 
cal orbitals. The retarded self-energy matrix £]'"(£') = 
Tip a'^p'" {E) is the sum of retarded self-energies intro- 
duced by the interaction with the leads which determine 
escape rates of spin-cr electron into the electrodes (p, a). 

For nonintcracting systems described by the Hamilto- 
nian ((21), the lesser self-energy is expressed in terms of 
%-(i?) as 



ll<{E)^Y.'r^E)t;{E). 



(8) 



P,(T 



The level broadening matrix in the rotating frame 



f;iE) = -2ImS;(i?) = -2ImS; [E + s—j, (9) 



hLLi\ 



is obtained from the usual self-energy matrices''^ 'Ep{E) 
of semi-infinite leads in the laboratory frame with their 
energy argument being shifted by sHuj/2 to take into ac- 
count the "bias voltage" in accord with Fig. [Ijb). The 
distribution function of electrons in the electrodes of the 
rotating frame DC circuit is given by 
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and charge 



^ ^^'^ 1 + exp[{E - Ef + shuj/2)/kTy 
where s = + for spin-| and s = — for spin- J,. 



III. LOCAL AND TOTAL CURRENTS IN 
TWO-TERMINAL TIIFM TI DEVICES 



(10) 



(5) 



bond currents flowing between nearest neighbor or next- 
nearest neighbor sites i and j if they are connected by 



The spatial imaging of local spin currents has played 
an important role in understanding how QSHE'^^ or 
mesoscopic SHE^ due to intrinsic SO couplings mani- 
fest in experimentally accessible multiterminal devices. 
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FIG. 3: (Color online) Spatial profile of nearest-neighbor and next-nearest-neighbor bond spin currents (darker arrow means 
larger current) pumped from the middle FM region of length dpM = 11 into the left and right 8-ZGNRs of length dcNR = 11 
with: (a) zero or (b) non-zero intrinsic SO coupling for which GNRs act as 2D TI. (c) Spatial profile of bond charge currents 
corresponding to (b). Right panels show bond current distribution over a selected transverse cross section (vertical twisty line). 
The FM magnetization is precessing with frequency hui/e « 8.3 fjY and cone angle 9 = 10°. 



In Fig. [3{a), we first establish a reference local-current- 
picture of pumping for ZGNR|FM|ZGNR device with no 
SO coupling in ZGNR regions. Its transport properties 
are governed^i by the edge-localized quantum states in- 
duced by the topology of the chosen zigzag edgesi^'^* The 
non-zero SO coupling in Eq. ([T]) converts the ZGNR re- 
gions into a 2D TI by opening energy gap and by forcing 
their strongly-localized states^i2i to become spin-filtered 
and to acquire a linear dispersion around kxd = ir cross- 
ing the band gap-'^^ (SO coupling additionally suppresses 
their unscreened Coulomb interaction^^). 

Using the picture of such helical edge states within the 
TI regions of the DC device in Fig.[T]^b), whose spin and 
chirality is illustrated in Fig. [3{b), we can follow possi- 
ble Feynman paths of electrons in Fig. [31 For example, 
a spin-| electron from (L,i) electrode at a higher elec- 
trochemical potential /i^ can only flow along the top left 
edge ^ then it precesses through the FM island (because 
it is not an eigenstate of term in J^rot) =^ enters with 
some probability into spin-l edge state on the bottom 
right edge =J> finally, it is collected by {R,1) electrode at 



a lower electrochemical potential /i^ < The spin-j. 
electron injected by (i?,i) lead would retrace the same 
path in the opposite direction on the way toward (L,t) 
electrode. 

According to these straightforward paths, one expects 
no chiral edge currents around bottom left and top right 
edges of the device. However, these currents do exist in 
Figs. [3{b) and Efc) due to more subtle backscattering 
effects at the TI|FM interface. They are the consequence 
of the paths, clearly visible in Fig. [3l where, e.g., spin- J, 
electron from (L,i) electrode is reflected and rotated^^ at 
the TI|FM interface (where FM region breaks the time- 
reversal invariance) to flow backward through spin-^ edge 
state =J> then it propagates along the left-NM-lead|TI in- 
terface =^ it flows along the bottom left edge, FM island, 
and bottom right edge to finally enter into (i?,t) elec- 
trode at electrochemical potential /x^. This process can 
also account for the difference between spin currents fiow- 
ing along the top and bottom edges of a single TI region, 
which yields non-zero total spin current in Fig. 2) 

Another set of transverse current paths, conspicuously 
visible in Figs, and [3Jc), emerges around TI|FM 
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FIG. 4: (Color online) 

at each transverse cross section of the two-terminal 
GNR|FM|GNR (730 = 0) and TI|FM|TI (730 / 0) devices 
obtained by summing all bond spin currents shown in right 
panels of Figs. [Sja) and |3jb) , respectively. Inset shows I^^ 
for a longer device of the same width — 8. The FM magne- 
tization is precessing with frequency hw/e ~ 8.3 and cone 
angle 61 = 10°. 
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FIG. 5: (Color online) The spin Hall conductance 
GsH = if'/V of a voltage biased eVi — —eV2 = eV/2 four- 
terminal bridge as a function of the size of a clean GNR cen- 
tral region playing the role of the QSH insulator of width 
Nz and length dcNR. The GNR is attached to four metallic 
electrodes modeled as "lattice-matched" square lattices whose 
contact to the honeycomb lattice is illustrated in Figs. [T]and 
[2] The transverse electrodes 3 and 4 cover the whole bot- 
tom or top zigzag edge of GNR and act as the voltage probes 

{V3 = V4=0^h=l4^ 0). 



interfaces. These interfacial spin and charge currents are 
able to penetrate slightly into the bulk of the TIs^^ in 
contrast to infinite homogeneous (i.e., all TI) systems 
where current flow is strictly confined to sample edgesi^ 
Figure m plots the total spin current along the device, 
which is obtained by summing local nearest-neighbor and 
next-nearest-neighbor bond currents at each transverse 
cross section in Fig. [31 Although non-zero locally, the to- 
tal charge current through any cross section in Fig.[3ljc), 
including the left and right electrodes, remains zero. 
Since typical spin-relaxation lengths are much longer 
than the length scale ~ hvp/ ^ over which pumping de- 
velops,— we assume that FM island is clean so that non- 
conserved spin currents emerge throughout its volume. 
The magnitude of spin current pumped into GNRs is set 
around the GNR|FM interface^i^ and it is enhanced by 
the presence of helical edge states (for 730/7 ^ 0.1), con- 
trary to naive expectation that FM|T1 interface would be 
less transparent for spin injection. 



IV. LOCAL AND TOTAL CURRENTS IN 
SIX-TERMINAL TI FM|TI DEVICES 

To obtain sharp conductance steps in quantum trans- 
port calculations of the integer QHE in realistic (e.g., 
consisting of ^ 10^ carbon atoms^) six-terminal meso- 
scopic all-graphene Hall bars requires to avoid reflec- 
tion at the lead| sample interface by using magnetic field 
both in the sample and in the semi-infinite leads, as 
well as by employing high quality contacts.'^^ Thus, it 



is somewhat surprising that perfectly quantized Gsh 
was obtained in Ref. |26| for a four-terminal graphene- 
based 2D TI device where SO coupling is present only 
in the sample and whose metallic electrodes, for a cho- 
sen square lattice model, have propagating modes that 
couple poorly to evanescent and propagating modes 
within GNRs-?2. In fact, the corresponding contacts 
metallic-electrode|zigzag-GNR in those devices act as ef- 
fective disorder by introducing mixing of transverse prop- 
agating modes defined by the semi-infinite ideal leads 3^ 

Figure [H extends findings of Ref. [2^, obtained for spe- 
cial graphene ribbon aspect ratios [N^ = 4n, c?gnr = 
8n-|- 1, n e N) and leads attached to a segment of the 
ribbon edge, to confirm that Gsh = 2 x e/47r can be 
obtained in any sufficiently^! wide and long graphene 
ribbon attached to four metallic electrodes. Here the 
square lattice leads cover the whole top or bottom lateral 
graphene edge, while the longitudinal leads are attached 
at armchair edges as shown in Fig. [21 The square lat- 
tice of the leads is selected to be "lattice-matched" (lat- 
tice constant is the same as carbon-carbon distance in 
GNR as illustrated by Figs. [Hand [J) in order to reduce 
the detrimental effects of the contacts (i.e., the trans- 
mission matrix of the two-terminal device NM|GNR|NM 
has smaller off-diagonal elements than when "lattice- 
unmatched" square lattice leads are employed to model 
metallic electrodes) 3. 

In order to generate non-zero total charge current re- 
sponse of the 2D TIs, we attach additional four electrodes 
(Fig. [S]) at the top and bottom edges of the device in 
Fig. [T^a) of length dcNR — 65|o?fm = 32|dGNR = 65 and 
width iVz = 32. The two longitudinal electrodes and ad- 
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FIG. 6: (Color online) (a) The spin Hall conductance Gsh = /V of a voltage biased Vi — V2 — V {Vp — 0, p—3-6) six- 
terminal dcNR ~ 65|dFM ~ 32|dGNR ~ 65 device of width — 32 whose static FM magnetization is coUinear with the z-axis 
orthogonal to the plane of the device. Insets label the six metallic electrodes modeled as "lattice-matched" square lattices 
shown in Fig. [T] Panel (b) plots the same quantity Gsh when static FM magnetization is tilted by an angle = 10° away from 
the jz-axis. For the same unbiased {Vp = 0, p=l-6) TI|FM|TI device, whose FM magnetization is precessing with frequency 
hu>/e « 8.3 /iV and cone angle 9 = 10°, panel (c) plots the longitudinal spin conductance Gsl = e/^'/ftw and longitudinal 
charge conductance Gl ~ el2/huj, while panel (d) shows the transverse charge conductance Gt = eli^jhio. In panel (d) we 
emphasize that condition = —le, as the signature of charge Hall current (i.e., any inverse SHE), is satisfied only within the 
SO energy gap Aso where GNR regions act as 2D TI. Note that the convention for the sign of the currents is Jp > (< 0) and 
Ip^ > (< 0) for particles flowing into (out of) lead p. 



ditional four transverse electrodes (covering 10 edge car- 
bon atoms) are modeled on the same "lattice-matched" 
square lattice. Such six-terminal bridge, when biased by 
the voltage difference Vi — V2 = V, exhibits quantized 
Gsh shown in Fig. EUa), despite the fact that it is not 
identical to standard homogeneous QSH bridges^^^. since 
it contains FM island in the middle breaking the con- 
tinuity of helical edge states. Furthermore, we confirm 
that quantization of Gsh is independent of the angle by 
which static FM magnetization is tilted from the z-axis, 
as illustrated by Fig.[S]Jb). 

Since Hlo < Ep, we can use f^{E)- f{E) w huiS{E- 
Ep) at low temperatures for the difference of the Fermi 
functions presenlji^ in Eq. ([3]). This "adiabatic approx- 
imation"— is analogous to linear response calculations 
for biased devices, allowing us to define the longitudinal 
spin conductance Gsl = e/^^'/Sa; and transverse charge 



conductance Gt — elu/hw for the device whose FM mag- 
netization is precessing with frequency lo. In the unbi- 
ased {Vp = 0) six-terminal TI|FM|TI device, Fig. [^c) 
shows that Gsl vanishes when Ep is within the SO gap, 
while non-zero charge Hall currents emerge in transverse 
electrodes so that Gt ^ in Fig. [n^d). Moreover, trans- 
verse charge currents obtained in the regime of the QSH 
insulator (marked by Aso 7^ gap in Fig. [6]) are the 
signature of the inverse QSHE since only in this range 
of Fermi energies we get 1^ — — /g and J3 = — Js char- 
acterizing the usual charge Hall effect. However, we find 
that Gt — el^/huj is not quantized. Note that the same 
conclusions arc reached when the electrodes are made of 
the same TI as in the central region of the QSH bridge. 

Figure [6ljc) also shows that total charge current 
pumped into the longitudinal leads is zero within the 
QSH insulator regime (Aso 7^ 0), which is quite differ- 
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ent from the conventional inverse SHE driven by spin 
pumping from precessing magnetization into a topologi- 
cally trivial SO-coupled systems, such as the multitermi- 
nal 2DEG with the Rashba SO coupling. In the latter 
case, AC charge currents (with small DC contribution 
vanishing as the 2DEG size increase) are pumped into 
both transverse and longitudinal electrodes. Their time- 
dependence originates from spin non-conservation in the 
presence of the Rashba SO coupling, unlike in our case 
where charge and spin-S'^ currents are time-independent 
in both rotating and laboratory frames since spin is 
conserved. 



V. DISCUSSION 

The presence of extended chiral edge states in QH 
and QSH systems makes the analysis of transport mea- 
surements based on the Landuaer-Biittiker multiprobe 
formulasi^ particularly simple assuming homogeneous 
bridge, such as the TI central region attached to elec- 
trodes made of the same TI material^"^'^ or 2DEG and 
graphene in high magnetic field attached to the same 
type of electrodes in the same magnetic field.~ In those 
cases, the extended edge states are perfectly matched 
across the whole device and one can simply draw^ii^ a 
picture of allowed trajectories (or "Feynman paths"), for 
whom the quantum edge states serve as guiding centers, 
to count their contribution (as ballistic one-dimensional 
conductors) to current in a selected electrode. 

Following Ref. [l|, we show such paths for the inverse 
QSHE in Fig. [71(a) based on an abstract scheme where 
pure spin current is injected through the longitudinal 
leads due to different electrochemical potentials for spin-f 
and spin-| states, fi\ — fjL\ = eV and ^\ — ^2 — —eV 
driving two counter-propagating fully spin-polarized lon- 
gitudinal charge currents. In this case, continuity of the 
edge states and absence of any reflection between the 
leads and the sample ensures that conductance associ- 
ated with transverse charge current is quantizedt^ How- 
ever, the setup is unrealistic from the experimental view- 
point because it requires separate control of electrochem- 
ical potentials for spin-| and spin-| carriers within the 
same electrode. 

On the other hand, the presence of FM region in our 
devices breaks the continuity of edge states, thereby mak- 
ing the network of Feynman paths much more compli- 
cated. The analogous issues have been explored in exper- 
iments^^ on QH bridges where the gate electrode covering 
small portion of the central region introduces backscat- 
tering between spatially separated edge states, thereby 
requiring complicated trajectory network22, or spatial 
profiles of local currents^^ to explain non-quantized fea- 
tures in the longitudinal resistance. Furthermore, in the 
case of our six-terminal TI|FM|TI bridge, the breaking 
of time-reversal invariance at the TI|FM interface by the 
nearby FM island also enables spin-dependent reflection 
where incoming electron from a helical edge state has its 
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FIG. 7: (Color online). Phenomenological analysis of the 
inverse QSHE via selected possible Feynman paths of spin- 
polarized injected and collected electrons for the case of: 
(a) four-terminal bridge of Ref. fH with homogeneous cen- 
tral 2D TI region (attached to TI electrodes) where pure 
spin current is injected from the longitudinal leads as two 
counter-propagating fully spin-polarized charge currents; (b) 
six-terminal bridge comprised of inhomogeneous TI|FM|TI 
central region whose FM magnetization is precessing at fre- 
quency Lu and six attached TI electrodes are made of the 
same material as the TI islands within the central region. 
The depicted Feynman paths in (b) are based on the twelve- 
terminal DC device in the rotating frame (with all of its elec- 
trodes accepting only one spin species and with any pair of 
electrodes of opposite spin being at different electrochemical 
potentials Hp — = ftw), and they are supported by the 
NEGF-computed exact spatial profiles of local charge cur- 
rents displayed in Fig. [S] The blue "L-shaped" connectors 
between two helical edge states of opposite spin and direction 
of propagation signify spin-dependent reflectionS^ at TI|FM 
interfaces. 



spin rotated to be injected in the counter-propagating 
helical edge state along the same edge of the sample. 

Nevertheless, some of the important paths can be ex- 
tracted from spatial profiles of local currents in Fig. [3] 
for the two-terminal TI|FM|TI device or Fig. [S] for the 
six-terminal TI|FM|TI device. This allows us to explain 
all of the key results on the inverse QSHE driven by spin 
pumping shown in Fig.[5]for total terminal currents. The 
possible electron paths are easier to draw for the six- 
terminal TI|FM|TI device whose electrodes are made of 
the same TI, and can be understood using the picture of 
twelve effectively half-metallic FM electrodes in the ro- 
tating frame which try to inject their fully spin-polarized 
electrons into chiral spin-filtered edge states moving in 
proper direction. 

For example, Fig. [Tj^b) shows spin-| electrons starting 
from lead (3, |) at a higher electrochemical potential /Ltg 
to enter the right moving spin-| helical state on the top 
edge. The electrons have probability to penetrate into 
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FIG. 8: (Color online) Spatial profile of charge currents in six-terminal devices whose central TI|FM|TI region with precessing 
FM magnetization is attached to: (a) six TI electrodes identical to central TI islands [as assumed also in Fig. Wi^)] ; (b) six NM 
electrodes (as in Fig. [Gjl modeled on the "lattice-matched" square lattice. The possible Feynman paths corresponding to current 
profiles in (a), as well as the labeling of the six terminals, is shown in Fig. Wih). The inset below each spatial profile depicts 
the corresponding total spin I^' = jf^" and charge currents I^' = Jij computed at each transverse cross section of 
the device as one moves from the left to the right longitudinal electrode. 



the FM island where they precess and continue to prop- 
agate down (e.g., through interfacial spin currents shown 
in Fig. [land Fig. [5] along the right FM|TI interface) to 
finally enter the bottom helical edge state as spin-t par- 
ticles which allows them to be collected by the electrode 



(6, t) at lower electrochemical potential /^g — Ms — 
At the same time, incoming electrons from lead (3, J,) 
can also be reflected at the TljFM interface where ac- 
companying spin rotation^^ makes it possible for them 
to propagate back into the lead (3,t) which they enter 
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as spin-t electron with some probability. Portions of this 
path are retraced by spin-| electron from lead (6,|) flow- 
ing toward lead (3,t), where reflection and spin rotation 
at the right TI|FM interface generate additional current 
of spin-t electrons into lead 6 to produce total charge 
current Iq > whose sign is compatible with the exact 
calculations shown in Fig. ^d). These paths also ex- 
plain the absence of quantization of Gt conductance in 
Fig. since transport which generates non-zero lon- 
gitudinal charge currents is not confined solely to the 
hehcal edge states as in Fig. El^a). 

While the selected paths in Fig. [7] are compatible 
with complete (numerically exact) spatial profiles of lo- 
cal charge currents shown in Fig. [51^a), the profiles also 
contain trivial strictly edge paths carrying non-zero local 
charge and spin currents between leads 4 and 2 due to 
ballistic transport of spin- J, electrons from electrode (4, |) 
toward (2, |) with higher state occupancy then for spin-l 
electrons propagating from (2, t) toward (4, 1) electrode. 
Since the same local currents exists along the opposite 
edge connecting electrodes 6 and 2, the total spin or 
charge current flowing into the longitudinal electrode 2 
is identically equal to zero [see inset below Fig. [SJa)]. 
The same conclusion holds by symmetry for leads 1, 3, 
and 5 on the opposite side of bridge. Thus, the absence 
of any Feynman paths that would allow electrons from 
(1,1) and (2,1) electrodes to reach some other (p,t) elec- 
trode via propagation through helical edge states, com- 
bined with trajectories within FM or along the TljFM 
interface, confirms one of the key results in Fig. ^c) — 
absence of longitudinal spin and charge currents in the 
spin-pumping-induced inverse QSHE. 

We also note that the above discussion becomes more 
complicated for the TI|FM|TI bridges with six NM elec- 
trodes that were actually employed in Fig. [51 This is ex- 
emplified by Fig. [HJb) where new paths have to be take 
into account that allow electrons to propagate along in- 
terfaces between TI regions and the attached NM elec- 
trodes while also penetrating^^ into the bulk of NM elec- 
trodes where there are no helical edge states. In the case 
of leads 1 and 2, local currents around TljNM interface 
penetrate only slightly into the bulk of the NM electrodes 
and then return toward TI to ensure zero total spin and 
charge currents in the longitudinal electrodes, as demon- 
strated by the inset below Fig. MJ^)- 



VI. CONCLUDING REMARKS 

Using the mapping of time-dependent spin pumping by 
precessing magnetization to a multiterminal DC device 
in the rotating frame, we describe injection of thus gener- 
ated local and total spin currents into helical edge states 
of a topological insulator in the absence of any externally 
applied bias voltage. In the regime where the voltage bi- 
ased six-terminal TI|FM|TI nanodevice exhibits QSHE, 
the corresponding unbiased device with precessing FM 
magnetization generates charge currents in the transverse 



electrodes (characterized by I4 = — Jg and I3 = —I5), 
while bringing pumped total spin and charge currents in 
longitudinal electrodes to zero. Although the transverse 
charge conductance of such inverse QSHE is not quan- 
tized, these two responses to pumping can be used to 
probe the TI phase via unambiguous electrical measure- 
ments. The absence of quantization of transverse charge 
conductances was explained as the consequence of elec- 
tron propagation paths being composed of both simple 
segments guided by chiral spin-filtered edge states of fi- 
nite extent and more complicated segments through the 
FM island or around the TljFM interface where spin- 
dependent reflection accompanied by spin rotation can 
take place. Our analysis based on imaging of local (i.e., 
on the scale of the lattice constant) spin and charge 
transport reveals interfacial currents around TljFM in- 
terfaces that can penetrate slightly into the bulk of the 
TI or interfacial currents around TIjNM-elcctrode con- 
tacts which break the continuity of helical edge states. 

Although a topologically trivial system with the intrin- 
sic SO couplings (such as the Rashba spin-split 2DEG^') 
would also generate non-quantized transverse charge cur- 
rents in response to pure spin currents pumped by pre- 
cessing magnetization)^ this setup also pumps longi- 
tudinal charge currents 1^ This is completely different 
from the behavior of the proposed TIjFMjTI multiter- 
minal device. In addition, for SO couplings which do not 
conserve Sz spin (such as the Rashba one), both trans- 
verse and longitudinal pumped charge currents are time- 
dependentf^S as opposed to DC transverse charge Hall 
currents generated by our device. 

During the preparation of this manuscript we became 
aware of the theoretical proposal^i for a charge pumping 
device designed to induce electrical response of a 2D TI 
where two FM islands, one with precessing and one with 
fixed magnetization, are deposited on the top and bottom 
edges, respectively, of a QSH insulator attached to two 
electrodes. While the physical motivation leading to this 
device — fractional charge response to magnetic domain 
wall acting as external physical field — is different from 
ours, its operation can be easily explained from our four- 
terminal DC device picture (Fig. [T]) of spin and charge 
pumping in the rotating frame. That is, the top FM is- 
land (see Fig. 3 in Ref. [Sj) injects pure spin current into 
upper helical edge states of the QSH insulator, which is 
then partially converted into charge current flowing along 
the edges as explained by spatial profiles of local charge 
currents in Fig.[3]Jc). The role of the bottom FM island in 
this setup is to block the local charge current, which oth- 
erwise flows in the opposite direction along the bottom 
edge in Fig. ^c) leading to Ii ~ I2 ^ in our device 
setup. Such blocking is based on a simple observation 
that FM region whose magnetization is coUinear with the 
magnetization of the effective half-metallic electrodes in 
the rotating frame does not permit any transport through 
it that connects these electrodes. However, it is unclear if 
such two-terminal device can display perfectly quantized 
longitudinal charge current / = ew/27r (or, equivalently. 
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quantized conductance Gl = e/2/Sa; = e^//i in our nota- 
tion) in a realistic setup, which was conjectured in Ref. lS 
via qualitative arguments without treating the effect of 
TI|FM interfaces on adiabatic charge pumping through 
either the scatteringii or NEGF— approaches. 
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